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ABSTRACT 

This work studies the appearance of a Haldane gap in quasi one-dimensional anti- 
ferromagnets in the long wavelength limit, via the nonlinear cr-model. The mapping 
from the three-dimensional, integer spin Heisenberg model to the nonlinear cr-model 
is explained, taking into account two antiferromagnetic couplings: one along the 
chain axis (J) and one along the perpendicular planes (J±) of a cubic lattice. An 
implicit equation for the Haldane gap is derived, as a function of temperature and 
coupling ratio J±/J. Solutions to these equations show the existence of a critical 
coupling ratio beyond which a gap exists only above a transition temperature Tn- 
The cut-off dependence of these results is discussed. 



I. INTRODUCTION 



It is by now well establishedmi that the excitation spectrum of the spin 
1 antiferromagnetic Heisenberg model has a mass gap (the Haldane gap) in 
one dimension, whereas it is gapless in higher dimensions!! Experimentally, 
integer spin chains are realized in quasi one-dimensional compounds in which 
the antiferromagnetic coupling J along the chain direction is much higher 
than the transverse. coupling Jj_. The Haldane gap can then be observed via 
neutron scatteringa or electron spin resonanceo. It is the ratio R = J±/J 
which determines the degree, of 'quasi one-dimensionality' of the material. 
For instance, R is estimatedo tOrbe ~ 0.02 in CsNiCbj, while it is certainly 
much lower (~ 0.0006) in NENP.H In CsNiCbj one-dimensional behavior (i.e. 
the existence of the gap) is oberved above a critical temperature (the Neel 
temperature TV) of about 5A". On the other hand, one-dimensional behavior 
is observed in NENP at temperatures as low as can be reached. The suggests 
the existence of a critical ratio R c below which the system is one-dimensional 
in character, whatever the temperature. 

This was argued in Ref. |?], wherein corrections to the spin-wave spectrum 
were claculated in perturbation theory for the anisotropic Heisenberg model 
in two dimensions. More recently, Azzouz and DoucotEa have performed a 



mean-field theory analysis on the Heisenberg mi 
H = J^Si •8 j + J±J2 S 

(ij) (Im) 

wherein the first sum runs over nearest neighboi 
the second sum runs over nearest neighbors in tl 
the chains. They found a critical ratio given, ir 
by R c ~ se~ 27rs , where s is the value of the inte, 

In this work we will perform a similar ana- 
dimensional non-linear cr-model, and at finite t 
used has been described in Ref. ||, wherein the 
the gap was calculated for the purely one-dimer 

The value we will find for R c is very sensitivi 
but if we fix the latter with the help of the nur 
the purely one-dimensional Haldane gap at zei 
critical ratio R c ~ 0.026. However, we believe 1 
regarded with great caution. We also illustrate t 
temperature, of the zero-temperature gap, and 1 
of the gap for various values of R. 



II. THE ANISOTROPIC NONLIN 

The aim of this section is to show how the Heis 
be replaced, in dimension 3 and in the continui 
nonlinear cr-model, with Lagrangian density 



C 



1 

~9 



1 



-(dttti) 2 — v(d z m) 2 — Rv(d x ii 



wherein 



= 2 JasVl + 2R 



2a2 
s 



and where a is the lattice spacing; the field m is 
constraint alone makes the theory non trivial. R 
(0) and (0) may proceed to the next section. 
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The mapping from the Heisenberg Hamiltonian. to the nonlinear cr-model 
has been done explicitly before in various wayscrclj, but mostly in dimension 
1. Here we wish to show explicitly how the mapping is done in dimension 3, 
with anisotropic couplings. We first need to write down the action for a single 
site, using spin coherent states. We define a unit vector n such that_S = sn. 
Then the kinetic term is well-known to be the Wess-Zumino action:E3 



(4) 



Swz = dt K = s dr dt n • (<9 T n x d t n) 



Here the time t runs from to some finite period T and r is an auxiliary 
coordinate introduced in order to parametrize, along with t, the spherical 
cap delimited by the curve n(t) from t = to t — T. For more details, the 
reader may consult Ref. [l^. For our purpose, it is sufficient to know that the 
variation of the action upon a small change Sn is 



SSwz = s / dt Sn ■ (n x dtn) 
Jo 



(5) 



Next, we add interactions between spins on a cubic lattice T', with nearest 
neighbor vectors (x, y,z). Anticipating short range antiferromagnetic order, 
we scale the unit cell by a factor of 2 in every direction, thus obtaining a 
bigger lattice T whose points are labeled r (unprimed) and whose unit cell 
contains 8 sites of V: 



r' = ax + (3y + 7Z 



(6) 



wherein each of a, ft and 7 runs from to 1. The unit cell of T contains 8 
spins which may be described by 8 differents vector fields, as follows: 



S(r') = s(-l) Q+/3+7 n(r') 

n(r') = m(r) +a^'(-l) M+J/?+ ' C7 l y fe(r) 



(7a) 
(7b) 



i,j,k 



Again each of 7 and k runs from to 1; we have introduced 7 
deviation fields 1^, the primed sum meaning that the term i = j = k = 
is omitted, the latter being rather represented by the slowly varying field m. 
We have included a factor of a in the definition of Ujk in order to stress that 
deviations from short range antiferromagnetic order are assumed to be small. 

Now we must express the kinetic term and the Heisenberg Hamiltonian in 
terms of these new fields, at lowest non-trivial order in a, and then integrate 
out the deviation fields Ujk to obtain an effective continuum action for m. 
Let us begin with the kinetic term K, which is obtained by summing over 
spins the action (Q). Assuming that n is slowly varying, we use Eq. (g) to 
write 



^ = S £]T(-ir+^n(r'Mn 

rer a, 13 

(r' = r + ax + (iy) 
where we defined the difference 

S z n(r') — n(r' + z) — n(i 
and similarly for S x and S y . The definitions (0) 

,5,n(r + ax + /?y) = -2a^(-l)" 



Applying the identity 



(-l) a ^ = 25, 



we may write the kinetic term as 

K = 2 d as^l m (r) • (m(r) x c 
rer 

where d is the dimension of space, which we ke- 
notation is three-dimensional, in order to make 1 
In the continuum, this becomes 



K = a 1 - d s J dr l m • (m x 



The Heisenberg Hamiltonian may be expresset 
for instance, 



V z = -s 2 J z n ( r ') ' n ( r ' 



r'er' 



Up to an irrelevant constant, this is equal to 

\s\j z ]T (^n(r')) 2 



r'er' 



We have similar expressions for V x ,y wherein J z 
and 5 X:V . We then use Eq. (|Io|), along with 

S z n(r + ax + fty + z) = 2a9 2 m(r) + 2a 'S 
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to find 



V z = 2 d a 2 s 2 J z ^2 i 2 ^2 l iji + 21 ooi ' d z m+ (<9 z m) 2 j (17) 

along with similar expressions for V x and V y . 

In the continuum, we may therefore write the following Lagrangian density: 



C = a 1 d slm • (m x d t m 

-d 

hj 



s2j A 2 Y, l *U + 21100 ' d x m+(d x m)'- 

i,j 

a~ d s 2 J y | 2 ^ I?y + 21 io • <9 a m + (<9„m) 5 
a- d s 2 Jz{zJ2 l<ji + 21ooi • 9 z m + (5 z m) 2 



(18) 



We now proceed to the functional integration of the fields ly*. Since the 
action is quadratic in these fields, this amounts to substituting in the La- 
grangian the expression of these fields obtained from the equations of motion. 
The latter are 



lm = 



1 



1 



-m x dtm 



(19a) 



4as J x + J y + J z 

lioo = -|d x m loio = -^dyin l 00 i = -\d z va (19b) 
ljjfc = otherwise (19°) 

Substitution of these equations into the above Lagrangian density yields ex- 
actly Eqs. @ and (g), with J z — J, J x = J y = RJ and d = 3. Note that in 
the isotropic case (J„ = J), the Lagrangian obtained in dimension d is 



with 



2a d ~ 1 

v = 2Jas\ /f d g = Vd 



(20) 



(21) 



The above derivation is valid for any dimension greater than 1. In dimen- 
sion 1 a complication occurs since there is only one deviation field 1, whose 
equation of motion is instead 



1 



m x d t m 



(22) 



Substitution into the Lagrangian yields, in adi 
topological Hopf term, responsible for the diff 
integer and half-integer spin chains. The above ( 
no such term arises in dimensions d > 1. 

So far we have eluded the question of constra: 
are related by 8 constraints coming from the rel 
small deviations (a — > 0) these constraints are 



1 



m ■ 1 



ijk 



the second of these equations is compatible wi 
(|H^), implying that the above substitution pre 
The constraint m 2 = 1 is part of the definition i 



III. DERIVATION OF THE GA! 

The aim of this section is to find, starting fr> 
(||), an equation governing the temperature an 
the Haldane gap. Before embarking upon calci 
the results, obtained in the approximation wher 
to the inverse lattice spacing a -1 ) is much lar; 
or the gap A. Let Ao be the Haldane gap at z 
dimensional limit (R — > 0), and let us introduc 
parameters 



t = 



T 
A^ 



r = 



Then the reduced gap S obeys the following equ 



. ln(<5 2 



-S 2 InS 2 - 



8i 2 




where we have defined the special function 



G(y) = \y 2 I dz ^z 2 - 1 [coth(yz/2) 



(K\ is the modified Bessel function of order 1). 1 
will be devoted to the proof of Eq. ( |25| ) . 

A large part of this section parallels the calci 
We start with the er-model Lagrangian expi 
variable ip = m/^/g: 
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C = -{(d ti p) 2 



(d zV ) 2 -R(d ± <p) 2 -<j(<p 2 -l/g)} 



(27) 



Here we have set the magnon speed v to unity; it can be restored at the 
end by dimensional analysis; we have also introduced a Lagrange multiplier 
a whose role is to enforce the constraint <p 2 = 1/g. 

The strategy used to compute the Haldane gap is to calculate at one loop 
the effective potential V(a) for the Lagrange multiplier. This amounts to 
integrating the field <p, assuming a to be constant. We then look for a 
minimum in V{a). If such a minimum exists for a ^ 0, then the position a 
of this minimum is the mass squared of the triplet field ip, namely the square 
of the gap. 

At one loop, the derivative V'(a) of the zero-temperature effective potential 
in the Euclidian formalism is 



1 



Pi 



(28) 



The gap equation is then simply V'(A 2 ) — 0, where A = yfa. The exten- 
sion of this formula to a finite temperature T is obtained by replacing the 
integration over po by a sum over Matsubara frequencies lu u — 2imT: After 
summing over Matsubara frequencies, the gap equation becomes (J3 = 1/T) 



i = 3 
9 



dp 1 



coth — 
2 



Rp 2 + A 2 



1>- 



(29) 



This expression is meaningless without the introduction of a momentum cut- 
off A, which should be proportional to the inverse lattice spacing a' 1 . We will 
evaluate this integral by restricting it to a cylinder of height 2A and radius A, 
parallel to the p z axis in momentum space. Defining the integration variable 
x = i?Pj, we have 



An 2 R 



, (Ik I dx , 1 coth-VPT^TA^ (30) 

3.9 Jo Jo Vfc 2 + x + A 2 2 V V ' 

The integral over x can be done analytically: 

— = T [ dk (in sinh ^ yjk 2 + A 2 + A 2 i? - In sinh ^ Vk 2 + A 2 1 (31) 
35 Jo I 2 2 J 

The remaining integral may also be done analytically, in the limit A > A, T: 
Let us introduce the notation 



F(a,/3,A)= [ dk lnsinh^Vfc 2 + a 2 
Jo 2 

After integrating by parts, this function becomes 



(32) 



F(a, (3, A) = A In sinh- ^ A 2 + a 2 

2 2 J 



k 2 



The second term may be separated into a contr: 



A k 2 dk 
Vk 2 



IA 2 + Ia 2 (l + In^ 



10 VK'tQ 

plus a finite temperature correction, finite as A 
rDO k 2 dk 



= (coth^x/^T 
a 2 I 2 v 



n 



dz \/ z 2 — 1 [coth( 



= 2a 2 G{(3a) 
Summing up all contributions and expanding ar 

F(a,0,A — > oo) = \(3K 2 - A In 2+ \fia? (1 - 



This result, inserted into Eq. (|3lj), yields at last 

^ 2R (A 2 + RA 2 )(i-m A2 + RA2 



3.9 



-8T 2 



G 



4A 2 
T 



This is essentially the gap equation (P5|), exc 
ables (|4[) have not been introduced. To this e 
dimensional case, i.e., let us take the limit R - 
that 



dG(y) 
dy 



-yH(y) 



H(y) 



We then obtain the one-dimensional equation 

3gA 2 4A 2 V 1 

Since H(oo) = 0, the zero-temperature gap is A 
terms of the reduced variables (0) the finite t 
dimension 1 is then 

±ln6 = H(5/t) 

This equation also appears in Ref. ^. Finally, ■ 
duced variables, the three-dimensional gap equal 
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IV. DISCUSSION 

In this section we illustrate the consequences of the gap equation (p5|). 
That implicit equation may easily be solved numerically, given the expression 
( ptf ) for the function G in terms of modified Bessel functions. Fig. 1 shows 
the t dependence of 8 for various values of the anisotropy parameter r. 

At zero temperature, the gap equation reduces to 



(5 2 + r) ln((5 2 + r) - S 2 In S 2 - r = 



(40) 



This equation has the solution S = 1 in the limit r — > 0, and has a nonzero 
solution for sufficiently small r. Fig. 2 illustrates the numerical solution of 6 
vs r. Taking the limit S — > we see that the critical value of the anisotropy 
parameter is r c = e. This in turn corresponds to a critical value R c of the 
coupling ratio R = Jj_ / J. 

Beyond this value of r, a gap will appear only above some critical temper- 
ature t c (the Neel temperature). If we take the limit 8 — ► in Eq. (p5|), we 
obtain 



r lnr 



■8ti 



G[f- )-**/6 



= 



(41) 



wherein we have substituted the special value G(0) = tt 2 /6. This equation 
may be solved numerically for t c as a function of r, and the result is illustrated 
on Fig. 3. 

Let us now discuss the relation existing between the critical value r c and 
the microscopic parameters J and R — J±/J. The coupling J resurfaces if 
the magnon speed v is restored in the explicit value of the gap by dimensional 
analysis. The expression for Ao, the zero-temperature gap in the R — > limit, 
is then 



A = 2vAcxp-(47r 2 /3gA 2 ) = 4JAasexp- [47r 2 s/6(aA) 2 ] 



(42) 



The ratio Aq/J can then be related to the cut-off prescription, i.e., to t 
product Aa. If we accept the ratio obtained from numerical simulation, 
namely Ao = 0.41J, we are led by the above equation to the prescription 
Aa = 2.1. Notice however that the ratio Aq/J is extremely dependent on 
the value of Aa, because of the exponential factor. 

After restoration of v, the expression for the anisotropy parameter r is 



r = ii?cxp- [8vr 2 s/6(aA) 2 ] 
The critical value r c = e then translates into a critial ratio 
R c = 4cxp [l - 87r 2 s/3(aA) 2 ] 



(43) 



(44) 



The choice Aa = 2.1 then leads to R c = .026. 
seems too high when compared to coupling ratk 
als; for instance, the ratio R has been estimatei 
in which a finite Neel temperature is observet 
other compound of the ABX3 type have a hexa 
It seems that trying to find the correct value ol 
off prescription is dangerous. However, definite 
from a continuum theory, up to an overall scale 
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FIG. 1. Reduced Haldane gap 5 as a function of 
various values of the anisotropy parameter r: 0.01, 



FIG. 2. Reduced Neel temperature t c as a functii 
r = RA 2 /A 2 . 
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FIG. 3. Magnitude of the reduced Haldane gap 8 = A/A at zero temperature 
as a function of the anisotropy parameter r = i?A 2 /A§. The critial value of r is e. 
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